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Nonlinear superposition of direct and inverse cascades in two-dimensional turbulence 

forced at large and small scales 

Massimo Cencini/ Paolo Muratore-Ginanneschi,^ and Angelo Vulpiani^ 
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We inquire about the properties of 2d Navier-Stokes turbulence simultaneously forced at small and 
large scales. The background motivation comes by observational results on atmospheric turbulence. 
We show that the velocity field is amenable to the sum of two auxiliary velocity fields forced at 
large and small scale and exhibiting a direct-enstrophy and an inverse-energy cascade, respectively. 
Remarkably, the two auxiliary fields reconcile universal properties of fluxes with positive statistical 
correlation in the inertial range. 
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Turbulence represents a tantalizing nonequilibrium 
system characterized by cascade processes which, as typ- 
ical in statistical physics, strongly depends on space di- 
mensionality. In d = 3, kinetic energy, injected at large 
scales, goes toward smaller ones with positive and con- 
stant flux (« energy-injection rate e), until is dissipated 
by molecular diffusion [l| . Between the injection and dis- 
sipative scales, the energy spectrum behaves as a power- 
law E{k) ~ e^/'^fc"^/? In d =2, ideal (inviscid and un- 
forced) fluids preserve both energy -< v^ > /2 and square 
vorticity (enstrophy) ^ w^ ^ /2 (w=V x v). On this ba- 
sis, Kraichnan [2| predicted that sustaining the flow at a 
single scale tj {^kj ), with energy (enstrophy) injection 
rate e {ij — ekj), generates a double cascade of enstro- 
phy downscale (< ^/) and of energy upscale {>if). He 
also predicted two power laws for the energy spectrum: 
E{k)Kjf/^k~^ (but for log-corrections [3|) in the direct 
enstrophy cascade range; E{k)K.€^/'^k^^^^ for the inverse 
energy cascade. The direct cascade, with a positive en- 
strophy flux («77), ends at the dissipative scale. Whilst, 
in an unbounded domain, the inverse cascade proceeds 
undisturbed, with a negative energy flux {K,—e), unless 
large-scale friction stops it at a scale ^€/ Ij]. For a recent 
numerical study of the dual cascade see [5| . 

In 3(i-layers, as the atmosphere, both 3d and 2d- 
phenomenology can be relevant depending on the aspect 
ratio, the injection and observation scales |6r^. Air- 
craft measurements [9|, [lOJ of atmospheric- winds revealed 
that horizontal energy spectra at the troposphere end (at 
K, lOKm altitude) display two power-laws: E{k) oc k^^/'^aX 
wave-numbers in the mesoscales (wlO — 500Km); E{k)oc 
fc-3 at synoptic scales (w500-3,000Km). Though, 2d 
phenomenology should dominate at scales larger than the 
troposphere thickness [Uj , measured spectra display the 
steeper and shallower power-laws in reverse order with 
respect to Kraichnan's scenario. To complicate the pic- 
ture, the energy flux seems to be positive at 10 — lOOKm 
12| . suggesting a Sd-like direct energy cascade, though 
the involved scales may be too large. 

In the 2d-framework, on which we focus here, sev- 



eral explanations for the observed spectra have been pro- 
posed. Interpreting the synoptic —3 spectrum as an en- 
strophy cascade, forced by instabilities of the horizontal 
motion at the planetary scale (^lO'^Km) [ll|, the —5/3 
mesoscales spectrum may result: from a 2d-inverse en- 
ergy cascade forced by convection driven by thermal gra- 
dients in the troposphe re Il3l4l5| : or, less likely |lClll6|. 
from gravity waves [lllll8|. Interpreting the —5/3 spec- 
trum as a convection-driven 2ci-inverse energy cascade, 
the steeper synoptical spectrum may result from large- 
scale coherent structures due to forcin g | l 9| or to energy 
condensation at the planetary scale [19|, l20j. Interest- 
ingly, such coherent motions may mask the inverse cas- 
cade inducing a positive energy-flux 19], which could 
explain observations [12|. Other proposed mechanisms 
(not discussed here) consider 3d or quasi-2(i scenarios ac- 
counting for stratification and other effects [21r.23j . 

In this Letter we focus on 2d turbulence forced at 
two, well separated, scales with the aim of understanding 
the interplay of oppositely directed cascades in the same 
range of scales. We thus consider the 2d incompressible 
Navier-Stokes equation, which for the vorticity uj reads 

dtLo+v-Vuj = -Upi-A)PLj~ag{-A)-iuj + f^+fi, (1) 

t)=V^V= (dyi/jj—dxip) is the velocity, and ip the stream 
function {A'ip — —uj). The hyperviscous (hypofriction) 
term removes enstrophy (energy) at small (large) scales 
generalizing standard dissipation p — 1 (Ekman friction 
(7 = 0). In direct numerical simulations (DNS), such gen- 
eralizations provide for extended inertial ranges. We use 
p = 8 and q^l; tests with different values have been also 
performed. The forcings /^ and fe act at separate scales 
L^i, injecting energy at scale i and enstrophy at scale 
L, at independent rates e and rj, respectively. We used 
two independent random, zero-mean Gaussian processes 
restricted to a narrow band in Fourier space k € [kn, ki2] 
(with i = i,L) centered at kz, ^ 1/L and k^ ~ l/£, with 
correlation ^fi{k,t)fj{k' , t)>-^F,5ij5{t-t')5{k+k')Q{k- 
kii)Q{ki2~k), Q being the Heaviside step function. In the 
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FIG. 1. (color online) Average energy spectrum Ev(k) (sym- 
bols). _Ei_2(fc) labels the energy spectra associated with the 
auxiliary velocities i;i,2 and suggests a direct (inverse) cascade 
for vi (1^2), see text. Arrows mark the forcing bands at /cl ~7 
and ki~240. Dotted lines show the slopes —3 and —5/3. In- 
set: energy 11" (fc) and enstrophy 11" (fc) fluxes normalized by 
their average values (e and rj, resp.). DNS of Eqs. (|H3|I were 
done with a standard 2/3-dealiased, pseudospectal method 
with 1024^^ collocation points. 
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FIG. 2. (color online) Dynamical fluxes of (a) energy no.^j(fe) 
and (b) enstrophy TlQ.jj{k) (defined in [231). Grey symbols in 
(b) refer to passive scalar flux 11®, see text (Cfr. Ea. (|10|l 'l. 
The solid curves show energy and enstrophy fluxes obtained 
by integrating ([2} with /l = (a) or /^ = (b), their su- 
perposition to IIo;oo and IIoioo, resp. further demonstrates 
the cascades overlap when both forcings are present. No- 
tice that no;22 ~ no;oo < and IIq;!! « making evident 
that V2 are the degrees of freedom associated with the inverse 
cascade. Similar considerations apply to vorticity, though 
IIo.ii > rEo-oo, which will be discussed later in relation to 
Fig.Bl 



sequel, we ignore the behavior at scales larger (smaller) 
than L (£) as damped by hypofriction (dissipation). 

In Rcf. 16], to explain the observed universality of 
atmospheric spectra, it was conjectured that the two 
sources may not be independent. However, large and 
small scale excitations originate from different physi- 
cal processes JJ| characterized by separate timescales 
(small-scale convection being faster) , bearing their inde- 
pendence. Within the independent-source model, spec- 
tral universality can be ascribed to that of the inverse- 
cascade yi. In IJ, |l6| it was also hypothesized that 



oppositely directed cascades could not coexist without a 
sink between the forcing scales. Lilly [15], using closure 
theories, showed that there is no need of such sink. Mal- 
trud and Vallis 2J| made, as far as we know, the unique 
numerical study of Eq. ([l}, providing evidence of two 
overlapping cascades. This is confirmed by Fig. [T] which 
shows that the energy spectrum E^(k) displays the ba- 
sic features of the atmospheric one. Moreover, enstrophy 
n"(A:) and energy 11" (fc) fluxes are constant with oppo- 
site signs (n"(fc) > and 11" (fc) < 0) meaning that direct 
enstrophy and inverse energy cascades superimpose, ap- 
parently undisturbed, in the same range [/c/, : kg]. We 
show below that, remarkably, the superposition is real- 
ized maintaining non-trivial correlations between the de- 
grees of freedom associated with the two cascades. 

To scrutinize this superposition we propose a decom- 
position able to disentangle the two cascades, by evolving 
in parallel with Eq. ([T]) two auxiliary equations 

dtu^, + V ■ Vu:= -^pi-Afuj,- a,i-A)-''uj,+ f, (2) 
dtoje + V ■ Vuji^ ~i^p{-A)PLUi - aq{-A)-''uJi+ fi , (3) 

where v is the same as in Eq. ([T]). Since /l and fe in 
(I2EI) are the same realizations of the forcings in ([1]), 
LOj^^e are two "active" pseudo-scalar fields [25|, such that 
Lo = uj^+LOe and v = v^+v^ [26|, with v^^i = V-^ip^^i 



(Aip^^i = —ujj^j) incompressible. In the sequel, we adopt 
the notation v = Vq, lu = u!o, v^^i = Vi^2 and ui^^i = u]i^2- 

The behavior of the energy spectra (Fig. [IJ and 
the fiuxes (Fig. [2|) associated with the auxiliary fields 
suggest the identification of the components i;i,2 as 
the carriers of the degrees of freedom mainly asso- 
ciated with, respectively, the direct and the inverse 
cascade. This observation can be substantiated by 
means of the Karman-Howarth-Monin (KHM) equation 
i {dt + a,(-A)-9)C,fe(r,i) + f,fe(r,t) - F,k{r,t) = 
^daSg jT {r, t) (as customary, we assumed translation, ro- 
tation and parity invariance). The KHM equation links 
the 3'"'^ order structure-tensor. 



'5"r"^ (r) =<5vr {r)Sv;^ {r)Sv^^ (r)^ 



(4) 



{5vi{r) = Vi{r,t) — Vi{0,t)), to the correlation functions 
of velocities {Cjk{r,t)—^Vj{r,t)-Vk{0,t))~) and forcings 
{Fjk{r) =-<fj{r, t)/fc(0, t)>-) and to the dissipative terms 
£jkir,t) = 2iyp -< V>f(r,i)V>f(0,t)^. Kraichnan's 
theory is equivalent to the following three hypotheses 
|28| : existence of steady state for Galilean invariant sta- 
tistical indicators; smoothness at finite viscosity; absence 
of velocity dissipative anomaly. Using these hypotheses, 
a careful analysis along the lines of 
in the range £ <S^ r ^ L, the expansion 



3 m Hal justifies 



S^r'ir)^En^osr'''P.^k{r'^^A 



{iljk]:n} 



5"."fc} 



■ P^Jk {A{,yk]:n}} r'^'r^^r'^" + o{r^) , 



(5) 



with shorthand notation for summation over cyclic per- 
mutations, Pijk{Oijk} = Oijk+Ojki+Okij for any O. 
Eq. ([5]) holds true strictly for i = in the ideal limit of 
infinite volume at vanishing hypofriction (aq — 0), with, 
furthermore, the explicit prediction for the coefficients: 

A{o[jk]:o}+Pajk {^{o[jfc]:o}}=2e((5jo+<5j2)(4o+42) (6) 
A{o[jk]:i} + iPojk {A{a[jk]:i}} ^ kv i^30+Sji){Sko+Ski)-{T) 
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structure functions: (a) for veloc- 



ity 'S'ooo('') compared with the prediction 



(sohd red hue) 
and Soi2{r); (b) for vorticity —Zo;ii{r) = ~r ■ Zo;ii{r)/r and 
Zo;i2{r), the black lines show linear behavior in r. The nega- 
tive linear behavior of Zo-ii{r) links to the direct cascade of 
^ Wi y; (c) Ci2{r) =^ LJi(0)cj2(j') >- vs r. The behaviors of 
Zo;i2{r) (which should vanish in the ideal limit) and of Ci2('') 
are evidences of the sensitivity of directly cascading degrees 
of freedom to the hypofriction, see text (Cfr. Eq. (O). 



Out of the ideal limit, the constants Auijk].n} = j^{i[kj]:n} 
depend on the full statistics of the solution of Eqs. ([I]- 
|3]). Moreover, for i — 1,2, no dynamical constraints can 
be imposed to fix the constants ^{i[jfc]:{o,i}}i however 
the expansion ([5|) can still be justified using parity in- 
variance and the incompressibility of the fields Vi. The 
quantities iSojfc are thus associated with dynamical fluxes, 
while Sijk and S2jk only provide statistical information, 
and will be dubbed "statistical" fluxes. Eqs. ([5][7]) recover 
the 3'' -order longitudinal structure function derived in 
Ref. gi|: 



'-'OOOV'') — '""l''a2''a3'-'000 v)/''' 



ler+ifjr^, (8) 



which is numerically very well satisfied (Fig. |3^). In- 
terestingly, other quantities deviate from the ideal limit 
prediction (e.g SQi2{r) = is not reproduced by data, see 
Fig. |3^). As discussed below, similar issues are present 
for 3'"'^-order quantities associated with uji {i — 0,1,2). 
The deviations from ideal-limit predictions should be as- 
cribed to the finiteness of the simulation domain, where 
the presence of hypofriction stymies the derivation of 
neat relations such as ([7]). The A^....iy's become coupled 
to the infra-red component of the kinetic energy hinting 
at a more non-local and less universal behavior of the en- 
strophy cascade. This fact relates to non-trivial correla- 
tions existing between the degrees of freedom associated 
with the two cascades. The importance of such corre- 
lations can be quantified by comparing the "statistical" 
fluxes n^:^ (fc) and 112.'^ (fc) (with j = 0, 1, 2), defined in 
[13, with the "dynamical" fluxes H^.'^^/s (for j = 0, 1, 2). 
As shown in Fig.|4l two phenomena stand out. 

Whilst njJ^oQ ~ 0^.22 < (Fig. Hi) validates the iden- 
tification of V2 as the degrees of freedom associated with 
the inverse cascade, the left panels of Fig. |4] show that 
ni;oo ^ (n5'.22 < 0) with intensity comparable to that 
of nj.Qo (n2.22) and thus indicates that Vi contributes to 
the inverse cascade of the total field Vq- This may appear, 
at first glance, surprising in consideration of the observed 
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FIG. 4. (color online) Fluxes of energy and enstrophy defined 
in [23]: (Left) energy fluxes lii-jj of Vj [j = 0, 1,2 from top 
to bottom) due to the transport by Vi (i = black semifiUed 
circles, i = 1 red empty circles, i = 2 blue filled circles); 
(Right) Enstrophy fluxes TLf.jj-. panels, symbols and colors 
follow the same convention of Left panel. 



absence of flux of kinetic energy of t>i (n^.^i ~Oforanyi). 
The non-intuitive and, possibly, non-universal behavior 
of the statistical energy fluxes H^ (« = 1, 2) is, however, 
a consequence of their dependence on the full statistics. 
Similar considerations apply to enstrophy fluxes, with the 
roles of Vi and V2 exchanged (Fig. fright). 

The second phenomenon appertains to the relative 
intensity of enstrophy fluxes. The ideal-limit energy- 
balance predictions ^JE^ suggest XIq.qq ~ Hg.^^ > with 



n; 



0;22 



0. The right panels of Fig. 2] confirm the latter 
prediction coherently with the interpretation of V2 as the 
carrier of the inverse cascade degrees of freedom. They 
also show that the "dynamical" enstrophy fiux XIo-n is 
significantly enhanced with respect to Hq.qq (see also 
Fig. EJj). The phenomenon can be regarded as a con- 
sequence of the alleged fact that the main contribution 
to Vi comes from the degrees of freedom of Vq undergoing 
the direct cascade and of the stronger sensitivity of this 
latter to non-local effects. In order to substantiate the 
claim we inspected the equations governing the correla- 
tion functions CfJr) =^uji{0)ujj{r)y. From the energy 
balance in the presence of hypofriction together with the 
assumption of inverse cascading V2 (i.e. ng.22 ^ with 
ag (-A)«C^2(0) = Fi =<ll>) we derived that 



no;ii - n;^;oo = l«. (-A)nQ - cri)(o) - f, 



0{kL) 



.gi-Ayc^.m + oikL) 



(9) 



for some 7 > 0, as enstrophy is bounded [28[. Eq. (|9]) 
evinces the non-local nature of the enstrophy transfer. Its 
positive value implies that the correlation €^2 {r) must be 
also positive in the inertial range, as confirmed in Fig. [3]:. 
In terms of the vorticity structure functions Zi-^jk{r) =< 
Svi{r)Sujj{r)Sujk{r)y, Eq. dH) translates into Zo-i2{r) = 
r ■ 2o-i2{r)/r > with linear dependence on r (as shown 
in Fig. 12)3), while it should be in the ideal limit. To 



further elucidate the enhancement of Hq.^^ with respect 
to Hq.oo, it is instructive to liken the ultra-violet behavior 
of the vorticity oji to that of a passive scalar 9 evolving 
in the same velocity field, as in Ref. J29J : 

dtO + v- Ve= -iyp{-A)P9- ag{-A)-'^9+ f (10) 
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Here, the large-scale forcing /^ has the same statistics 
of that used in ([2]) but with independent realizations, 
making 9 a passive scalar 25|. Flow incompressibility 
brings about a scalar energy -<9'^>-l2 cascade to smaller 
scales, and similarly to ([S]), by energy balance, we have 



n'^-n: 



0:00 



«,(-5.')''(Co"o-0(0)/2 + 0(fcLr (11) 



for some positive 7', where C^(r) =^ 9{r , t)9{0, t) y . We 
expect the difference in ([TT|) to be positive as the infra-red 
content of Cqq is "fattened" by the inverse energy transfer 
a.t k < k]^. As shown in Fig.[2j), the comparison between 
the passive scalar energy flux, H^, with the enstrophy 
flux, IIq.qq validates this prediction. The inference is that 
the same phenomenon underlies the positive sign of (J9l). 

Summarizing, we showed that 2d-turbulence sustained 
by a large and a small scale forcing gives rise in between 
the sources to an inertial range where a direct and an 
inverse cascade co-exist and overlap. We also showed 
that there exists a natural decomposition of the Navier- 
Stokes field compartmentalizing the degrees of freedom 
associated with the direct and inverse cascade in two aux- 
iliary velocity fields, obtained considering a single large 
and small scale source, respectively. Although these aux- 
iliary fields satisfy energy balance relations as if they 
were independent, actually they are not and exhibit non- 
trivial correlations pinpointed by the inspection of third 
order statistics of "statistical" objects, evading the en- 
ergy balance relations. In contrast to the settings used 
here, realistic forcings in the atmosphere should be time- 
correlated, with the large-scale excitation being slower 
than the small-scale one. Provided the forcings are inde- 
pendent with separate spatial and temporal scale the pic- 
ture presented here should remain essentially unaltered. 
However, a slow large-scale forcing may induce coherent 
motions that, as argued in Ref. [19], can change the sign 
of the (total) energy flux and thus mask the inverse cas- 
cade process, which may explain observations 12[. 

We conclude with a theoretical remark. The decompo- 
sition in terms of auxiliary velocity fields here proposed 
can be straightforwardly generalized to 2d Navier-Stokes 
equations with an energy input distributed over differ- 
ent scales. In this perspective, the cascade overlap of the 
two sources model, here investigated, evinces the physical 
mechanism for why Kraichnan theory applies also in the 
presence of power-law sources and, consequently, for the 
inability of renormalization group approach to correctly 
predict Navier-Stokes energy spectra |30| . even in what 
may seem a priori a perturbative regime. 
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